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Abstract 

We propose a new parametrization of the deceleration parameter to study its time-variation 
behavior. The advantage of parameterizing the deceleration parameter is that we do not need to 
assume any underlying theory of gravity. By fitting the model to the 157 gold sample supernova 
| la data, we find strong evidence that the Universe is currently accelerating and it accelerated in 

the past. By fitting the model to the 115 nearby and Supernova Legacy Survey supernova la data, 

o . 

, the evidence that the Universe is currently accelerating is weak, although there is still a strong 

°. 

evidence that the Universe once accelerated in the past. The results obtained from the 157 gold 



sample supernova la data and those from the 115 supernova la data are not directly comparable 
because the two different data sets measure the luminosity distance up to different redshifts. 

We then use the Friedmann equation and a dark energy parametrization to discuss the same 
problem. When we fit the model to the supernova la data alone, we find weak evidence that the 
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Universe is accelerating and the current matter density is higher than that measured from other 
experiments. After we add the Sloan Digital Sky Survey data to constrain the dark energy model, 
we find that the behavior of the deceleration parameter is almost the same as that obtained from 
parameterizing the deceleration parameter. 
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I. INTRODUCTION 



Astronomical observations suggest the existence of dark energy which has negative pres- 
sure and contributes about 2/3 of the critical density to the total matter density of the 
Universe within the framework of Einstein's general relativity . While a wide va- 

riety of dynamical dark energy models were proposed in the literature Q, there are also 
model-independent studies on the nature of dark energy by using the observational data. In 
particular, one usually parameterizes dark energy density or the equation of state parameter 
w(z) of dark energy F° r the parametrization of w(z), we need to determine 

Q m o in addition to the parameters in w(z). Furthermore, one needs to assume the validity 
of general relativity in all these studies. 

Since supernova (SN) la data measures the luminosity distance redshift relationship, it 
provides a purely kinematic record of the expansion history of the Universe. It is possible 
to probe the evolution of the Hubble parameter or the deceleration parameter by using SN 
la data without assuming the nature and evolution of the dark energy jlll . Irj . llS^ . For 
instance, by parameterizing the deceleration parameter q(z) — q + qiz, Riess et al k| 
were able to study the kinematics of the universe. However, it was soon realized that such 
a parametrization cannot re-produce the behavior of the cosmological constant Q|. An 
alternative parametrization is a piecewise constant acceleration with two distinct epochs 



ll| , for which it was found that SN la data favors recent acceleration and past deceleration. 
Recently Shapiro and Turner (ST) applied several simple parametrization of the deceleration 
parameter q(z) to the 157 SN la data j3] to study the property of q(z) 3]. They found 
that there is little or no evidence that the Universe is presently accelerating, and that there 
is very strong evidence that the Universe once accelerated. The advantage of parameterizing 
q{z) is that the conclusion does not depend on any particular gravitational theory. The 
disadvantage is that it will not give much direct information on the cause of an accelerated 
Universe. 

The conclusion arrived in was based on a simple three epoch model of q(z), in which 
the function q(z) is not smooth. Since the current SN la data is still sparse, the division 
of the data to three different redshift bins may not be a good representation of the data. 
Therefore, the conclusion based on this technique needs to be further studied. Following 
ST, we propose a simple smooth function of q(z) which we believe is more realistic and then 
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apply the observational data to get the behavior of the deceleration parameter. After we 
are sure that the Universe experienced acceleration, we assume that the acceleration is due 
to the presence of dark energy and use a simple dark energy parametrization to study the 
property of dark energy in the framework of general relativity. 

Specifically, the paper is organized as follows. In section II, we study the property of 
q(z) by fitting the parametrization q(z) = 1/2 + (qiz + <72)/(l + z) 2 to the 157 SN la data 
and the 115 nearby SN la and the Supernova Legacy Survey (SNLS) SN la data compiled 
in jl6]. In section III, we apply the parametrization w(z) = w + W\zj (1 + z) 2 to study the 
property of q(z) and the nature of dark energy. The Sloan Digital Sky Survey (SDSS) 3] 
and the Wilkinson Microwave Anisotropy Probe (WMAP) data [^| are also combined with 
the SN la data in our analysis. In section IV, we conclude the paper with some discussion. 



II. THE CURRENT ACCELERATION 



From the definition of the Hubble constant H(t) = a/a and the deceleration parameter 
q{t) = —a/(aH 2 ), we have 



H(z) = H exp 



[1 + q(u)]dln(l +u) 



(1) 



where the subscript means the current value of the variable. So if we are given a function of 
q(z), then we can find the evolution of our Universe without applying any particular theory 
of gravity. ST considered a simple three-epoch model [3] 



g for z < z t , 

Q{ z ) = { qi for z t < z < z e = 

q2 = 0.5 for z > z e = 0.3. 
By using Eq. (JTJ), for the simple three-epoch model we get 

H (l + z) l+ i» 



0.3, 



H(z) 



H (l + z t y^ (ig 



1+91 



1+91 



l+Z 



1+92 



for z < Zt, 
for z t < z < z e 
for z > z e . 



The parameters in the model are determined by minimizing 

2 [Vob S (Zi) - fi(Zi)} 2 
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where the extinction-corrected distance modulus n(z) = 5 log 10 [dz,(z)/Mpc] + 25, is the 
total uncertainty in the SN la data, and the luminosity distance is 

ST found that a long epoch of deceleration is consistent with the 157 gold sample SN la data 
at the 10% level, and they concluded that there is little or no evidence that the Universe is 
presently accelerating. Although this model is very simple, the functions q(z) and H(z) take 
different forms at different epoches. Especially, q is not continuous. Because the current 
SN la data is still sparse and the decomposition to only three redshift bins is not a good 
representation, the conclusion derived from the model may not be robust. In this paper, 
we would like to use both the 157 gold sample SN la data and the 115 nearby SN la and 
the SNLS SN la data compiled in |lf|. In this 115 data set, there is no SN la with redshift 
0.101 < z < 0.249 which is around the redshift z t in Eq. (J2J). Therefore it is not a good 
idea to fit the simple three-epoch model to the 115 SN la data, so we proposed a simple 
two-parameter function 



= n + H , .A 2 ' ( 6 ) 



1 qiz + q 2 

2 (1 + z) 

to fit the 115 SN la data. Note that q(z) —>■ 1/2 when z ^> 1 and go = 1/2 + q 2 , so the 
parameter q 2 gives the value of go- The behavior of q(z) in this parametrization is quite 
general except that q(z) —>■ 1/2 when z ^> 1 which is consistent with observations. If q± > 
and q 2 > 0, then there is no acceleration at all. It is also possible that the Universe has been 
accelerating since some time in the past. If q\ < and q 2 > —1/2, then it is possible that 
the Universe is decelerating and has past acceleration and deceleration. In other words, this 
model may have the same behavior as the simple three-epoch model. The values of gi and 
q 2 and the behavior of q(z) can be obtained by fitting the model to the observational data. 
Substitute Eq. ® into Eq. JTJ), we get 



H(z) = H (l + zfl 2 



exp 



<Z2 q\z 2 - q 2 



(7) 



2 ' 2{l + z) 2 _ 

Fitting the model to the 115 SN la data, we get x 2 — 113.65, gi = — O.8I2J an d = 
— 1.1 5I035, here the given error is the la error. By using the best fitting results, we plot 
the evolution of q(z) in Fig. ^ From Fig. ^ we see that go < 0, i.e., the Universe is 
currently accelerating, contains over 96% of the probability. At the 3cr level, it is possible 
that g > 0. So it is possible that the Universe is decelerating now although the evidence is 
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FIG. 1: The evolution of q{z) = 1/2 + (q 1 z + q 2 )/{l + zf by fitting it to the 115 SN la data. The 
solid line is drawn by using the best fit parameters. The dotted lines show the \a error, the dashed 
lines show the 2a error and the dotted dash lines show the 3<7 error. 

not strong. It seems that we have strong evidence that the Universe had acceleration in the 
recent past. We see that the transition redshift when the Universe underwent the transition 
from deceleration to acceleration is z t = 0.95^11 at the lcr confidence level. However, the 
redshift range for past deceleration is around 0.2 which is in the same range that there is no 
SN la data. To make a more solid conclusion, we also fit the model to the 157 gold sample 
SN la data and found that \ 2 = 174.07, q 1 = 1M±%™ and 92 = -l.59t.oil- We plot the 
evolution of q(z) by using the fitting results in Fig. |21 

From Fig. |21 we see that q(z) < for < z < 0.2 contains over 99.7% of the probability. 
This result suggests that there are strong evidence that the Universe is currently accelerating 
and the Universe once accelerated. Fig. El shows that the transition redshift Zt = 0.36^Qog. 
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FIG. 2: The evolution of q{z) = 1/2 + (q x z + q 2 )/(l + -z) 2 by fitting it to the 157 gold sample SN 
la data. The solid line is drawn by using the best fit parameters. The dotted lines show the la 
error, the dashed lines show the 2a error and the dotted dash lines show the 3a error. 

III. DARK ENERGY PARAMETRIZATION 



Although the above simple model fits the data well and it tells us that the Universe 
experienced acceleration, but it hardly tells us anything about the property of dark energy. 
In order to make connection with dark energy, we apply Einstein's general relativity and 
work the problem in the usual way. We work on the dark energy parametrization 

W\Z 

w{z) =w + ( 1 + ^ 2 > (8) 

because Wq may be positive and w(z) —> wq when z ^> 1. This parametrization may give a 
currently decelerating universe. It is valuable to mention that the above parametrization is 



just a Taylor expansion in the scale 
wi(a/ao) 2 , so it has limitations too 



'actor a(t) to the second order, w(z) = w +Wx(a/a ) — 
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The dimensionless dark energy density is 



Sl DE (z) = n DE0 (l + zf( 1+Wo) exp [3^72(1 + zf] . (9) 
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FIG. 3: The evolution of q(z) by fitting the model w(z) = wo + w\z/(l + z) 2 to the combined 157 
gold sample SN la and SDSS data. The solid line is drawn by using the best fit parameters. The 
dotted lines show the la error, the dashed lines show the 2a error and the dotted dash lines show 
the 3a error. 

Therefore, we get 

Ono + (1 - Ono)(l + z) 3m [l + 3w + 3w lZ /(l + z) 2 } exp [3 Wl ,z 2 /2 (1 + z) 2 } 
q[Z) 2[fi m0 + (l-fi m0 )(l + ^exp[3^ 2 /2(l + z) 2 ]] ' 1 j 

So g = (^mO + (l-^mo)(l + 3u>o))/2. When w = 0, we get q = 0.5 > 0. When w = -0.4, 

we can get q > for reasonable value of fl m o- The model was used to fit the combined 157 

gold sample SN la , the SDSS and the WMAP data in [10]. Unfortunately, there was an 

sign mistake in w a in the results. We fit the model again to the combined 157 gold sample 
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0.04) w 



-1.72; 



-0.72 
-0.66 



SN la and SDSS data and we find that x 2 = 172.84, fi m0 = 0.26 
and w\ = 6.23^5;72. Combining the contour data of and Wi with Eq. (JTUJ), we plot the 
evolution of q(z) in Fig. |H1 From Fig. |2l it is evident that q(z) < for < z < 0.2 contains 
over 99.7% of the probability. Fig. H shows that z t = 0.30l{j;gjj. 



Now let us fit the model to the 115 SN la data alone, we get the best fitting values 



r 



113.36, Q m0 = 0.3918;^, w = -0.971^ and w 1 = -4.4l^ 5 . So w = is within the 



la error of the best fitting result, which means that for the fitting to the 115 SN la data 



TABLE I: Summary of the best fit parameters to the combined 115 SN la and SDSS data. 





w 


Wi 


qo 


x 2 


0.27 ±0.04 


-0.99±g;|| 


-0.30±1;^ 


-0.58 


113.55 


0.28 ±0.03 


-0.4 


-6.2±i;° 


0.068 


117.13 


0.29 ±0.03 





-lOAtH 


0.5 


122.80 



alone, it is possible that the Universe is experiencing deceleration now. If we fit the model to 
the 157 gold sample SN la data alone, the result remains the same. In fact, if we set wo = 
and fit the 115 SN la data, we get the best fitting values x 2 — 115.53, f2 m o = 0.461q;q9 and 
Wi = —21.9^5 9 . The result is consistent with the result of setting w free. Note that the 
best fitting value of f2 m0 is higher than that obtained from other observations. To fix this 
problem, we add the SDSS data to the SN la data. So we add the term (A - 0.469) 2 /0.017 2 
to Eq. (J3J). The parameter A = 0.469 ± 0.017 measured from the SDSS data is defined as 

1/3 

(11) 



_ H \/Q m0 



0.35 



#(0.35) \J H(z 



0.35 / f - 35 dz 



The best fitting results to the combined 115 SN la and SDSS data for the model (JHJ) are 
X 2 = 113.55, tt m0 = 0.27 ± 0.04, w = -0.99^;g and w x = -0.30l|;^. Combining the 
contour data of wq and W\ with Eq. (jlOj) . we plot the evolution of q(z) in Fig. 0] From 
Fig. HI we see that the evidence that go < is weak and there is strong evidence that the 
Universe once accelerated. Fig. HI shows that z t = 0. 771^36 • Of course, we may add the shift 
parameter derived from WMAP data to fit the model. By adding WMAP data, we found 
that the conclusion remains the same. We also find that wq = —0.4 is within 2a error and 
Wo = is within 3a error. For u7o = —0.4, we find the best fitting results are x 2 = 117.13, 
Qmo = 0.28 ± 0.03 and w\ = — 6.2t\'®. For u> = 0, we find the best fitting results are 
X 2 = 122.80, Q m0 = 0.29 ± 0.03 and W\ = — 10.4j^g. These results are also summarized in 
table 1. 

In Fig. El we plot the relative magnitude for the three cases discussed above with respect 
to the ACDM model with Q m0 = 0.27. The higher the redshift, the bigger the difference. 
While it is easier to distinguish the current deceleration models with the ACDM model with 
more accurate observational data, it is more difficult to distinguish the general dynamical 
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FIG. 4: The evolution of q(z) by fitting the model w(z) = wo + w\zj(l + z) 2 to the combined 115 
SN la and SDSS data. The solid line is drawn by using the best fit parameters. The dotted lines 
show the la error, the dashed lines show the 2a error and the dotted dash lines show the 3<7 error. 

dark energy model with the ACDM model. 
IV. DISCUSSION 

By fitting the parametrization q(z) = 1/2 + (q x z + g2)/(l + z ) 2 to the 157 gold sample 
SN la data, we find that q(z) < for < z ^ 0.2 within the uncertainty of 3a level. Recall 
that q(z) < means that the Universe is in the acceleration phase at redshift z. We also 
get z t = 0.95lQ5g at the la level. If we fit the parametrization to the 115 nearby and SNLS 
SN la data, we find that qo > at the 3a level and strong evidence that the Universe once 
accelerated. We also find that z t = 0.36^q q§ at the la level. Note that these conclusions do 
not depend on any particular theory of gravity. 

To study the role that general relativity might play, in this paper we also considered the 
Friedmann equation together with the dark energy parametrization w(z) = wq+W\Z / (1+z) 2 . 
After applying them to the 157 gold sample SN la data or the 115 SN la data, we find that 
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FIG. 5: The dependence of relative magnitude-redshift relation upon parameters wq and w%, 
relative to the ACDM model with f2 m o = 0.27. All parameters are the best fit values to the 
combined 115 SN and SDSS data. The bottom panel plots the curve wq = —0.99 and w\ = —0.299 
again so that the scale can be seen more clearly. 

go > at the la level and the best fitting Q m0 ~ 0.4, which is higher than that determined 
by other observations. In the parametrization q(z), we have two parameters q± and q2- 
When we use the dark energy parametrization, there are three parameters Q m o, wo and w%. 
With the addition of one more parameter apart from the assumption that the acceleration 
is due to dark energy in the framework of general relativity, we expect that the constraints 
will be loose. Therefore the dark energy parametrization is then fitted to the combined SN 
la and SDSS data because the combined data will break the degeneracies among the three 
parameters. If we fit the dark energy model to the combined 157 gold sample SN la and 
SDSS data, we find that q(z) < for < z ^ 0.2 contains over 99.7% of the probability. 
We also get z t = 0.30±g;§| at the la level. If we fit the dark energy model to the combined 
115 SN la and SDSS data, we find that go > at the 2a level and strong evidence that the 
Universe once accelerated. The transition redshift is found to be z t = 0.77j^ ; 3 g at the la 
level. Whether we use the parametrization q(z) or the dark energy parametrization w(z) 
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to fit the SN la data, the evidence that the the Universe once accelerated is very strong. 
These results confirm the conclusion that there is strong evidence that the Universe once 
accented obtained in Q. Ur resnlts also sn gg est that „ > 0.2. Althongh the upper 
bound on z t is too large by using the parametrization q(z), z t is more tightly constrained to 
be z t < 1.0 by using the dark energy parametrization w(z). 

The situation about go is more subtle. It depends on the model and the data. If we fit 
the models to the 157 gold sample SN la data, we find strong evidence for go < 0. But if 
the models are fitted to the 115 SN la data, the conclusion is different. The discrepancy is 
caused by the difference of the redshift range the data probed. The 157 gold sample SN la 
data probes much deeper in the redshift (z ~ 1.7) than the 115 SNLS data which stops at 
z ~ 1. For the dark energy model, the evidence for current deceleration looks promising if 
the model is fitted to the SN la data alone. The fitting also gives larger value for Q m o- When 
the SN la data is combined with the SDSS data for the dark energy model, we get reasonable 
value for Q m o and the sign of go is uncertain. The uncertainty of the question whether the 
Universe in accelerating now is directly related to the value of wq- If wq —0.6, then the 
Universe is accelerating currently. If Wo < —1, then the dark energy behaves like a phantom 
and it is not explained by the quintessence field. In fact, different models mean different 
external priors, so it shows that whether the Universe is currently accelerating depends on 
external prior. 

In conclusion, we confirm that there is strong evidence that the Universe once accelerated. 
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